Time-reversal symmetry relation for nonequilibrium flows 
ruled by the fluctuating Boltzmann equation 
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A time-reversal symmetry relation is established for out-of-equilibrium dilute or rarefied gases 
described by the fluctuating Boltzmann equation. The relation is obtained from the associated 
coarse-grained master equation ruling the random numbers of particles in cells of given position and 
velocity in the single-particle phase space. The symmetry relation concerns the fluctuating particle 
and energy currents of the gas flowing between reservoirs or thermalizing surfaces at given particle 
densities or temperatures. 
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O ■ I. INTRODUCTION 

a 

Since 1872, Boltzmann's equation has provided the main paradigm of our understanding of irreversible phenom- 
ena In isolated systems such as a dilute gas in a container, the iJ-theorem established by Boltzmann with his 
equation shows that any velocity distribution for the particles irreversibly converges at long time towards the Maxwell 
velocity distribution characterizing the thermodynamic equilibrium Q ■ This relaxation towards equilibrium is gener- 
ated by successive binary collisions between the particles composing the gas. They are described as in chemical kinetics 
by the mass action law, requiring that the rate of binary collisions be proportional to the concentrations of particles 
of the corresponding velocities in the volume element where the collisions happen. For this reason, Boltzmann's 
equation is nonlinear as the mean-field kinetic equations describing macroscopic chemical reactions, except that the 
particle velocities are not macroscopic observables 0. Nevertheless, finer observables closer to the microscopic level 
O ' of description and, especially, fluctuations remain outside the framework of Boltzmann's theory. 

Since the forties, a fluctuating Boltzmann equation has been proposed which rules the local velocity distribution 
I ■ function as a random variable, much in analogy with Langevin's stochastic equation (34l3| . This formulation provides 
, a description closer to the microscopic level and, thus, more suitable to understand the properties of the fluctuations. 
^\ • Although they are known to be time-reversal symmetric at equilibrium because of the principle of detailed balancing, 
C*~) ' few results are available about fluctuations in nonequilibrium steady states. 

The purpose of the present paper is to establish a time-reversal symmetry relation valid out of equilibrium for the 
' fluctuating Boltzmann equation. Such relations have been the subject of the so-called fluctuation theorems, which 
have been obtained in particular for stochastic processes such as diffusion processes ruled by Langevin's equations and 
their associate Fokker-Planck equation or continuous-time Mar kovian jump processes p^ - [l7| . In gen eral, a fluctuation 
theorem holds for all the currents flowing across an open system in nonequilibrium steady states [18| . Such a theorem 
implies the second law of thermodynamics and it allows us to deduce generalizations of the Green-Kubo formulae 
and Onsager reciprocity relations for the nonlinear response properties [ljl [2(|. Our aim is here to extend these 
fundamental results to the dilute and rarefied gases ruled by Boltzmann's equation (2l| - [23j . We consider such gases 
flowing under nonequilibrium conditions in pores, pipes, or other ducts between several reservoirs, as illustrated in 
Fig.ffl 

The paper is organized as follows. In section|TTl Boltzmann's equation is introduced as a mean-field equation in open 
geometries with gas-surface interactions. The symmetries of the collision kernel and, in particular, its time-reversal 
symmetry are discussed. In section IIII1 the one-particle phase space is partitioned into cells and the coarse-grained 
master equation is obtained for the probability that the cells contained certain particle numbers. In section IIV1 the 
diffusive approximation of the coarse-grained master equation is shown to lead to the fluctuating Boltzmann equation 
including the terms due to the contacts with the reservoirs. The time-reversal symmetry relation is proved in section 
fVl where the corresponding fluctuation theorem is established. Conclusions are drawn in section IvTl 



(a) 



(b) 




FIG. 1: Schematic representation of dilute or rarefied gases flowing through pores or pipes between reservoirs at given temper- 
atures and particle densities. In gas mixtures, there are several particle species i and so many particle densities. 



II. BOLTZM ANN'S EQUATION AS A MEAN-FIELD KINETIC EQUATION 

At the microscopic level of description, the motion of the particles composing the gas is ruled by Newton's equations 
for their positions and velocities {r n (t), v n (i)}. The particle distribution function is defined as the density to find one 
particle with the position r and the velocity v at the current time t: 

/(r, V) i)5ji[r- r«(i)] S [v - v n (t)] (1) 

n 

for monoatomic particles without internal rotation or vibration. Given that the initial conditions of the particles are 
distributed according to some probability distribution for the whole system, the particle distribution function Q is a 
random variable which may fluctuate. We may also consider its average value over the given probability distribution: 

</(r, v, t)> ee J2(S [r - r n (t)] 5 [v - v n (t)]> (2) 

n 

which is expected to smoothly vary in space and time. 



A. In the bulk of the flow 



For a general classical system, the time evolution of the one-particle average distribution function ([2]) is ruled 
by the first equation of the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy [2J-127J. The other equations of this 
hierarchy are ruling the many-particle average distribution functions. In the dilute-gas limit, the many-particle 
distribution functions factorize in terms of the one-particle function. The time evolution of the one-particle function 
is determined by the binary collisions and their differential cross section, which is obtained using classical scattering 
theory. Supposing that the particles are incoming every binary collision without statistical correlation as stated in the 
famous "Stosszahlansatz" , Boltzmann's equation is deduced, which is a closed equation for the one-particle average 
distribution function of the following form: 



8(f) | v 8(f) { F 8(f) 
dt dr m <9v 



= J rfv 2 rfvi dvX Vl ,v 2 |vi,v^) ((/() (}> 2 ) - (fx) </ 2 » (3) 



where {/) = (/i) is the average distribution function at the position r = ri and the velocity v — vi of the first particle 
involved in the binary collision j23 - [27| . F = F(r) is an external force field, which includes the repulsive forces of the 
walls of the duct. In general, the transition rates have the following symmetries: 

time-reversal symmetry: w (vi, v 2 |v' 1 , v 2 ) = w (— v^, — v 2 | — vi, — v 2 ) (4) 

space-orthogonal symmetry: w (vi, v 2 |v' 1 , v 2 ) = w (R • Vi, R • v 2 |R • v' l7 R • v 2 ) (5) 
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where R is a matrix belonging to the orthogonal group 0(3) including spatial rotations and reflections (2(| [27j- We 
also have the 

inversion symmetry: w (vi, V2|v^, v 2 ) = w (v^v^lvi, V2) (6) 

which is the symmetry under the time-reversal transformation combined with the spatial inversion of orthogonal 
matrix R = — I where I is the identity matrix. 

The transition rates are related to the differential cross section Cdifi of the binary collisions as 

w (vi, v 2 |vi, v 2 ) = CTdiff 8 (vi + v 2 - v[ - v 2 ) S (v 2 + v 2 - vf - v 2 2 ) (7) 

where the delta's express the conservation of linear momentum and kinetic energy in every binary collision. Setting 
V = (v[ + v 2 )/2 and u' = v' x — v 2 and integrating over the velocities dv^dv'2 = dV'du' = dV'u' 2 du'dfl, Boltzmann's 
equation ([3]) takes its more usual form: 

<>.j) v _ ay) + f _ ay) = r ^ dn ^ ^ ^ ((/ , } (/ , } _ (A) {h)) (g) 



dt dr m d 



v 



in terms of the differential cross section and the solid angle element dfl. For a gas of hard spheres of radius a, this 
cross section is given by Odiff = a 2 - 

In the absence of force field F = 0, Boltzmann's equation is known to admit stationary solutions given by the 
overall equilibrium Maxwell-Boltzmann distributions: 

m/3\ 3/2 ^ 1 _ 2 



/eq(v) = y-^-J 71 ex P("^ e ) with e = 2 TOV ( 9 ) 

for particles of mass m at the uniform temperature T = (/cb/3) -1 and density n. Since the transition rates preserve 
the conservations of linear momentum and kinetic energy in the binary collisions, the time-reversal symmetry Q also 
implies the property of reciprocity: 

w(vi,V 2 |vi,V 2 ) /oq(vi) / cq (v 2 ) = W(-Vi,-V 2 | - Vl,-V 2 ) /cq(vi) / cq (v 2 ) (10) 

We notice that the same property still holds in an inertial frame moving at the velocity Vo where the equilibrium 
distribution is given by Eq. ([9]) with v replaced by v — Vo. Moreover, the property of reciprocity as well as the 
stationarity of Boltzmann's equation for the Maxwell-Boltzmann distribution hold locally at every space point in 
the system. These two features constitute the basis of what is called the local thermodynamic equilibrium in this 
system [24j |. 

In the following, the form §3§ of Boltzmann's equation will be used because its symmetries ^ and §5§ are more 
explicit. 

B. At the boundaries 

For open systems such as those depicted in Fig. [1] there are several types of boundaries. 

1. At the reservoirs 

First, there are the inlets from the reservoirs, by which gas is also evacuated in the opposite directions. The contacts 
with the reservoirs contribute to Eq. by terms of the form: 



d(f) 



Of 



J2 I d 2 S S(r - r s ) (n • v) 9(n • v) [/ s (v) - </(r, v; t))] (11) 



where 8(x) is Heaviside's function, r s is the position of a point on the surface s separating the system from the 
corresponding reservoir, n is a unit vector normal to the surface s and directed inside the system and 

fm(3 s \ 3/2 1 
is( v ) = -7T 1 ™s exp(-^ s e) with e = -mv 2 (12) 



2tt / n r ' 2 
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is the Maxwell-Boltzmann distribution function of the gas incoming through the surface s from the reservoir at the 
temperature T s — (fce/Ss)" 1 and the particle density n s . The contacts with the reservoirs are carried out at surfaces 
which are assumed to be far enough in the pipes for the flow profile to be independent of the details of the boundary 
conditions. 

These inlets drive the system out of equilibrium if the reservoirs are at different temperatures or particle densities. 
For instance, in the case of the system of Fig. [TJi with two reservoirs s — L, R and one species of particles, there are 
two possible differences and thus two thermodynamic forces also called affinities which are the control parameters of 
the nonequilibrium constraints. These affinities are defined as 

thermal affinity: A E = /3 R — /3 L (13) 

n L /3 3/2 

chemical affinity: An — /?lA*l — /?rMR — m |tj (14) 

Ma 

where /xl,r are the chemical potentials of the monoatomic gas in the left- and right-hand reservoirs. We notice 
that more than two affinities are defined for an open system in contact with more than two reservoirs at different 
temperatures or particle densities as illustrated in Fig.QJ}. If all the affinities are vanishing, the system is at equilibrium 
and Boltzmann's equation admits the corresponding Maxwell-Boltzmann distribution as stationary solution. 



2. At the surfaces 



Now, we consider the interactions of the gas with solid surfaces, on which the particles are scattered. In the simplest 
model, the gas particles undergo elastic collisions on the surface, so that the outgoing velocity is given by 

v' = v-2(n-v)n (15) 

where n is a unit vector normal to the surface at the point of collision. Accordingly, the velocity component perpen- 
dicular to the surface changes its sign: n • v' = — n • v. Moreover, the kinetic energy is conserved because v' 2 = v 2 . 
Therefore, the contact of the gas with the surface does not change its temperature or density, as if the gas was in a 
force field F(r) that is very steep at the surface. In this case, the contribution of the surface to Eq. ^ would be of 
the form 



d(f) 



dt 



d 2 S S(r - r s ) (n • v) 0(n • v) {{f[r, v - 2(n ■ v) n; t]) - </(r, v; t))} (16) 



In general, the interaction of the gas with the surface of a solid depends on many different aspects. The scattering 
process of particles with surfaces have been much studied and many processes are known besides elastic collisions: 
adsorption, desorption, transport on the surface, transport into or from the bulk of the solid, or possible reactions f28j . 
The solid surface is typically at some temperature that may differ from the temperature of the gas and a transfer 
of energy can happen at the surface or inside the solid. Here, the solid forming the container is supposed to have 
a high enough thermal conductivity so that its temperature is uniform. Under certain thermodynamic conditions 
of pressure and temperature, gas particles may diffuse into the solid and form a stable thermodynamic phase such 
as solid hydrides, oxides, or nanoporous composites. Under nonequilibrium conditions, the solid could thus become 
a sink or a source of gas particles, in which case the difference of chemical potentials of the particles between the 
gas and the solid would add an extra affinity. For simplicity, this situation is not considered here and only thermal 
exchanges between the gas and the solid are envisaged. Surface adsorption and desorption may also occur but the 
surface is supposed to be in local thermodynamic equilibrium with the gas so that the thermodynamic state of the 
surface (in particular, the particle coverages) is thus stationary. Under such circumstances, the contribution of the 
surface to Eq. ^ would be of the form 



d(f) 



dt 



s 



= / d 2 S S(r - r s ) [ f dv'|n-v'|p r (v|v') (/(r, v'; i)> - (n • v) 0(n ■ v) </(r,v;t)) 

Js Unv'<0 



(17) 



where p r (v|v') is the probability density that a particle impinging the surface at the position r with the velocity v' 
and n • v' < will be scattered to the velocity v such that n • v > [23j. We notice that the probability density 
does not depend on the position r if the surface is homogeneous and the scattering process is the same everywhere. 
In general, this function is normalized according to 

p r (v|v') dv = 1 if n-v'<0 (18) 

nv>0 
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and it satisfies the following two properties. The hrst one guarantees the preservation of equilibrium at the temperature 
T s of the surface: 

[n-v|/ s (v)=/ Pr(v|vO |n-v'|/ s (v')dv' (19) 

Jnv'<0 

where / s (v) is the Maxwell-Boltzmann equilibrium distribution (fl2"|) at the temperature of the wall. The second is 
the property of reciprocity: 

|n • v'| / s (v') p r (v|v') = |n- v| / s (v) p r (-v'| - v) (20) 



which is implied by the time-reversal symmetry of the underlying microscopic dynamics and the condition of local 
thermodynamic equilibrium of the surface at the temperature T s of the Maxwell-Boltzmann equilibrium distribu- 
tion / s (v) HI. 

The special case (fT5|) of elastic collision is recovered for 

p r (v|v') = 5 [v' - v + 2(n • v) n] (21) 

uniformly on the whole surface, which satisfies all the aforementioned properties. In general, energy exchange happens 
during the lapse of time between the adsorption of the gas particles and their desorption. If a local thermodynamic 
equilibrium at the temperature of the solid is assumed, the particles are expected to be desorbed with the corresponding 
Maxwell-Boltzmann distribution. A very simple model achieving this condition is provided by the kernel 

p r (v|v') = |n • v| exp(-/3 s e) with e = ^mv 2 (22) 

which is independent of the other velocity v' and which also satisfies the required properties. More realistic gas-surface 
kernels have been considered and discussed in the literature [23| . In principle, the form of the probability density 
should be determined by studying the quantum scattering process of the particles with the surface on the basis of 



their atomic structure [28 1. 



Remark. There exist further possible situations if surfaces at in relative motion one with respect to another, as 
in Couette- Taylor flows. In such cases, the relative velocities of the surfaces are extra parameters controlling the 
nonequilibrium driving of the flow. In some other circumstances, the linear momentum of the flow can be conserved, 
e.g., if the gas moves between parallel surfaces on which the particles undergo elastic collisions. In this case, the 
velocity distributions of the gas incoming from different reservoirs could have non-zero average velocities and an 
extra affinity is associated with the difference between the average velocities (29|. However, linear momentum is not 
conserved for typical geometries of the surfaces and gas-surface interactions so that the chemical affinity may be 
sufficient to characterize the nonequilibrium constraint on particle transport. 



III. THE COARSE-GRAINED MASTER EQUATION 



At the mesoscopic level of description, the quantities of interest are typically fluctuating. This is for instance the 
case for transport processes or chemical reactions where the macroscopic equations of time evolution are replaced at 
the mesoscale by stochastic equations of Langevin type and their associated Fokker-Planck equation S 12l.JT3l ■ 

Similar considerations have been applied to the Boltzmann equation and related kinetic equations [3l4l3l| . The idea 
consists in supposing that the fluctuating distribution function ([T]) is ruled by a stochastic process. In order to obtain 
the master equation of this process, the one-particle phase space of the position and velocity variables (r, v) can 



be partitioned into cells centered around the phase-space points (r a ,v a ). 
number of particles in the cell a at the current time t is defined by 



These cells are of volume Ar 3 Av 3 . The 



N a (t) 



f /(r,v,t) drdw ~ /(r a ,v a ,i) Ar 3 Av 3 

J a 



(23) 



Since the distribution function / is fluctuating, these particle numbers are random variables undergoing jumps every 
time a particle exits or enters the corresponding cell. 

The probability to find {N a } particles in the cells {a} at the current time t is denoted as 



P = P(N) = P({N a }) 



(24) 
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When transitions happen, particles are exchanged between the cells so that it is useful to introduce the rising and 
lowering operators 

S± 1 $(..., JV a , ...) = #(..., JV a ±l,...) (25) 

acting on any function $(N) = <f>({iV a }) by adding or removing one particle in the cell a. These operators can be 
written as 

E^ = exp(±4r) (26) 



. ~dN a 

which implies that they are adjoint to each other: 

(S± 1 ) t =^ 1 (27) 
The master equation ruling the probability distribution (l24l) takes the following form 0, [HI EH 

£ Wl™ (E^E? ~ l) N P P + £ W% pa {E-^E^EfE? l) JV^P (28) 



dt 



Xp X\Xp(7 



with the transition rates due to the free flow (F), the contacts with the reservoirs (R), and the interactions with the 
surfaces (S): 



<f S = Wl p + W& + Wl (29) 



and those due to the binary collisions W 



Xppa " 

The time evolution of the average particle numbers 



(N a ) =Y,N a P(N) (30) 



N 

and their higher statistical moments is ruled by a hierarchy of equations that can be deduced from the master 
equation (l28l) [12j . For dilute gases, the second moments can be supposed to be factorized in terms of the averages 
themselves, in which case these latter quantities obey the following closed mean-field equations: 

!<*«>=£ (N p ) - (N a ) < Q RS + 2 E W ^p° W W 2 W E (N„) (31) 

p A p-pa \pa 

Although this mean-field kinetic equation is nonlinear because of the binary collisions, the master equation is linear 
and can be written as 

f = lp < 32 > 

in terms of the linear operator L given by Eq. (|2"5)l . We notice that this operator is composed of two terms, L = 
iv _l_ ^ one corresponding to the motion of independent particles and the other to the binary collisions. 

At equilibrium, the stationary solutions of the coarse-grained master equation (|28[) are given by multiple Poisson 
distributions of the following form: 

^^^^'^T with (^)=/ cq (v tt )Ar 3 A« 3 (33) 

a 

being the average values of the particle numbers in every cell a given in terms of the overall equilibrium Maxwell- 
Boltzmann distribution §§§ at the uniform temperature T — (&b/3) -1 an d density n. This multiple Poisson distribution 
represents the equilibrium state because it corresponds to the Maxwell-Boltzmann distribution ^ that is the equi- 
librium stationary solution of Boltzmann's equation ([3]). 
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A. In the bulk of the flow 

The transition rate of the events during which one particle is flowing in free flight p — > A from the cell p to the cell 
A is given by 

W lp ~ ^ ( U *P ' V ) 9 ( n *P • V ) ^ A ,r p + Arn Ap <5v A ,v p (34) 

where n\ p is a unit vector directed from the center r p of the cell p to the center r>, of the cell A [13, [Hj]. The 
Kronecker delta £>r p .r A -Arn Ap expresses the fact that the transition occurs between next-neighboring cells such that 
ta = r p + Ar n\ p . The other Kronecker delta 5 VpjVA means that the velocity remains constant during the free flight. In 
the case where the particle undergoes an acceleration, <5 Vp .v A should be replaced by the appropriate expression taking 
into account the force field F(r) acting on the particle. Since the one-particle motion is time-reversal symmetric, the 
corresponding rates have the symmetry: 

time- reversal symmetry: W Xp — WJt x t (35) 

where a T = (r Q ,,v a ) T = (r a ,— v Q ) is the time-reversed cell corresponding to the cell a. The symmetry (|35[) holds 
because n\ p = —n p \. Moreover, the rates ([34]) satisfy the property of reciprocity: 

Wl p (N p ) cq = AT (N x ) cq (36) 

where {N a ) eq denotes the average number of particles in the cell a at the equilibrium state (|33p . 
On the other hand, the rate of the transition pa — > Xp due to binary collisions is given by 



W c 



^— - [ dv 1 { dw 2 { dv\ ( dv' 2 ^(vi.valv'^v^) 5 Tli>rx 6 Tpirx 5 Ttr>Tx (37) 

J X J n Jp Jo 



These transitions do not change the positions but modify the velocities according to the collision rule. We notice that 
the symmetries (U)-© of the transition rates of Boltzmann's equation imply similar symmetries for the rates (|37[) : 

time-reversal symmetry: W^ pa . = W^ (T t a t a1 t (38) 

inversion symmetry: W^ prT = W° aXli (39) 

Moreover, since the binary-collision terms in the master equation (|28[) are summed over the four indices, the following 
symmetries also hold Q 

W^ pa = = Wg ap = WX P (40) 

The rates (|37)l also satisfy the property of reciprocity: 

WaU (^>«i ( n °^ = Kw* 1 ( N ^« ( N ^« ( 41 ) 

so that the coarse-grained master equation (|28l) preserves the local thermodynamic equilibrium. The multiple Poisson 
distribution Q33p is symmetric under time reversal because the Maxwell-Boltzmann distribution is invariant under 
velocity reversal v — > — v so that (N a ) = {N a T) at equilibrium. Therefore, the principle of detailed balancing 
holds, according to which the probabilities of every transition are equal to those of the opposite transition at the 
thermodynamic equilibrium [l2l [ill. [34|. In the limit Ar 3 Au 3 — > where the volume of the cells vanishes, we recover 
Boltzmann's equation © in the bulk of the flow from Eq. (l3"Tj) with the transition rates (|3"4l and ([3"7| . 



B. At the boundaries 



1. At the reservoirs 



For open systems, the partition into cells {a} should be carried out in a finite volume between the reservoirs. The 
particles are entering and leaving the system from and to the reservoirs. These specific motions are ruled by the flow 
operator i FRS . The key point is that this operator depends on the boundary conditions imposed to the system by 
the temperatures and particle densities of the flows incoming from the reservoirs. We first consider the simplest case 
of free flight between two reservoirs, as schematically depicted in Fig. [5J 
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FIG. 2: Schematic phase portrait of the free motion between two reservoirs from which the particles are incoming with different 
temperatures and densities. 



The time-reversal transformation (r, v) — > (r, — v) maps the trajectories with positive velocity coming from the 
left-hand reservoir onto those with negative velocity coming from the right-hand reservoir. The phase space between 
both reservoirs is partitioned into cells. The cells at the boundaries exchange particles with the reservoirs. The 
particles are dropped into the corresponding reservoir where they exit the open system. However, they enter into the 
system with a Maxwell-Boltzmann distribution at the temperature and density of the reservoir from which they come. 
Accordingly, the rate of the transition N — > N + 1 in the cell at the boundary of an incoming flow is determined 
in terms of the temperature and density of the nearby reservoir in analogy with the master-equation formulation of 
diffusion [32j. In the simple geometry of Fig. [2] the part of the flow operator for particles in free flight at the velocities 
±v acts on functions $({JVi,±}f = i) where iVj ± is the number of particles in the cell at the position i and velocity ±i>. 
The free flight from one reservoir to the other is partitioned into I cells of size Ar so that the cell index i runs from 
i = 1 to i = I. For the particles of velocity ±i>, the flow operator with the inlets from the reservoirs would be: 



L FR $ = 



v 

Ar 



i-i 



(=i 

i-i 



(E+X - l)JVi,_* + ^2(E+X lt _E7X - l)N i+h _<f> + [EJX - 1)(JV_) R $ 



(42) 



where (-/V + )l = JiAr 3 Av 3 is the average number of particles incoming with velocity +v from the left-hand reservoir 
at the temperature Tl and density while (JV_)r = f^Ar 3 Av 3 is the average number of particles incoming with 
velocity — v from the right-hand reservoir at the temperature 7r and density These average numbers can be 
expressed in terms of the equilibrium Maxwell-Boltzmann distribution function (|12[) in the corresponding reservoir 
s = L, R. For simplicity, only the particles moving at the kinetic energy e = mv 2 /2 in one direction are here considered. 

The average particle numbers at the reservoirs are related to the thermal and chemical affinities (fT5|) and (fT4]) 
according to 

A = eA E +A N = ln^±=\n^ (43) 
V*-)R JR 

The stationary solution of the master equation (1251 for non-vanishing values of the affinities gives the invariant 
probability distribution of the corresponding nonequilibrium steady state. 

At equilibrium, the stationary solution of the master equation (1251 is the multiple Poisson distribution ([33)) so that 
the temperatures and chemical potentials are equal and all the affinities are vanishing. 



2. At the surfaces 



In the case of elastic collisions on a surface, the particle follows a deterministic orbit which is piecewise straight. 
Let us consider an orbit coming from the right-hand reservoir, colliding on the surface and going back to the same 
reservoir. The speed ||v|| is conserved during the whole flight although the velocity direction undergoes the specular 
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reflection given by Eq. (fTS"]) . This is the case for every orbit of this kind. Accordingly, the motion along such orbits 
can be described as before by a partition into I cells of size Ar so that the cell index i runs from i = 1 to i = I. For 
an orbit interrupted by a specular reflection, the flow operator would be: 



Ar 



v 



1-1 



1-1 



(44) 



Comparing with Eq. (|4^|) . we see that all the terms are now of propagative type except the terms at the right-hand 
reservoirs, but they are the same as before. This means that the operator is of the same type as for a free flight from 
the right-hand reservoir back to itself and no difference in temperature or particle density manifests itself in this case. 

The situation is different in the case where gas-surface interactions induce energy exchange described for instance 
by Eq. (TTTj) with the simple thermal kernel (p?!?|) . To fix the ideas, the surface is taken as z — so that the normal 
unit vector is n = (0,0, f). ff we discretize position and velocity, the corresponding operator in the coarse-grained 
master equation is given by 



J2 S z , £ E A ^ PM V ') 4-v 1 - l) N rv $ (45) 



v'<0 



expressing the adsorption of a particle of velocity v' with v' z < and its desorption at the velocity v with v z > 
distributed according to a kernel satisfying the three properties (TTS)) . (fT9"|) , and (f!?U|) inducing a thermalization at the 
temperature T s of the surface. This operator has to be added to the free- flight operator L F . 

IV. THE FLUCTUATING BOLTZMANN EQUATION 

In the diffusive limit, a functional master equation of Fokker-Planck type is deduced from the coarse-grained master 
equation (|28j) . The diffusive limit is reached by considering cells containing sufficiently large numbers of particles that 
they can be described by Gaussian distributions while taking the limit of arbitrarily small cells to reach the continuum 
limit. The resulting functional equation of Fokker-Planck type rules the time evolution for the probability density 
functional that the random distribution function ([lj takes some given form at the current time t: V[f] — V[f(r,v)]. 
In the basic situation without contact with reservoirs or thermalizing surfaces, the functional master equation writes: 

§- t m = -jdrd, ^ A[f] V[f] + lj d r dv * dV , /(r;V /; /(r , v/) B[f] V[f] (46) 

with the deterministic kinetic rate 

Af] = -v • % - - ■ ^ + /dv 2 dvi dV 2 w (vx, v 2 |vi, v0 (f[ f 2 - h h) (47) 
or m ov J 

and the diffusivity coefficient 

B[f] = S(t-t') J dv 2 c?v 3 dv 4 w (vi, v 2 |v 3 , v 4 ) (/1/2 + /3/4) 

x [8(vi - v') + <5(v 2 - v') - <5(v 3 - v') - 5(vi - v')] (48) 

and where S/Sf denotes a functional derivative [13, EH, E3|- We notice that this functional master equation is of 
diffusion type as for the usual Fokker-Planck equations associated to Langevin stochastic differential equations. The 
average distribution function @ obeys an equation at the top of a hierarchy of moments equations. If the noise 
amplitude is low enough, this hierarchy may be truncated and Boltzmann's equation is thus recovered. To go 
beyond this mean-field approximation and consider the effects of the fluctuations, the functional master equation (|46[) 
should be used. 

If the system is open, the injection of particles from the reservoirs is described at the mean-field level by the term 
(fTTjl in the Boltzmann equation ©. For the coarse-grained master equation ([25)1 . this is taken into account by the 



10 



boundary terms in the flow operator (|42|) . In the diffusive limit, the following extra terms must be added to the 
deterministic kinetic rate (|4"T|) and the diffusivity coefficient (|4"5)l : 

a Js 

B R [f] = S(r - r') <5(v - v') £ f d 2 S S(r - r s ) (n • v) 0(n • v) [/.(v) + /(r, v)] (50) 

Similarly, the presence of a thermalizing surface also contributes by corresponding terms .4 s [/] and B s [/] . 

The analogue of the Langevin equation associated with the functional Fokker-Planck equation (|4"6")l is the following 
stochastic integro-differential equation of Ito type for the random distribution function JTJ) : 

+ I dv 2 dv[dv , 2 w(v 1 ,v 2 \v' i y 2 ) (f{f 2 -hf 2 ) + g(r,v,t) (51) 

ot or m ov J 

where g(r, v,i) is a Gaussian white noise satisfying: 

(ff(r,v,i))=0 (52) 

(g(r,v,t)g(T'y,t'))=5(t-t')5(r-r') J dw 2 dv 3 dw 4 w (v x , v 2 |v 3 , v 4 ) (/1/2+/3/4) 

x [5(vi - v') + «5(v 2 - v') - 6{v 3 - v') - «5(v 4 - v')] (53) 

as shown explicitly in Ref. [10]. In an open system, the contact with the reservoirs and thermalizing surfaces adds 
the contributions 

(, 9 (r, v, t) g(r>, v', t'))*' S = 8{t - t') B R > s [f] (54) 

to the variance (1531) of the noise. 

Since the functional master equation (|46l) can be viewed as the continuum limit of the coarse-grained master 
equation (|28[) , we shall use this latter which is more convenient to handle mathematically. 

V. THE TIME-REVERSAL SYMMETRY RELATION AND ITS CONSEQUENCES 

A. The theorem 

To fix the framework, we consider the open system with two reservoirs sketched in Fig. [T^, in a stationary state of 
affinities A = {Ae, An} given by Eqs. (ITBl and (fl4"l) . The transport of energy and particles between the reservoirs is 
characterized in terms of the instantaneous currents: 

+00 

jE{t) = <rke k S(t-t k ) (55) 

k— — 00 
+00 

3N(t) = <rkS(t-t k ) (56) 

k— — 00 

where {t^t^oo ^ s the sequence of random times when a particle enters or exits the open system at the left-hand 
reservoir with its energy e^. The sign of the exchange is taken into account by setting a = +1 if the particle enters 
and a = — 1 if the particle exits. 

We introduce the cumulant generating function of the fluctuating currents as 



Qa(A) = lim — In ( exp 

t— »oo t v 



-A - [ j(r)dr 
Jo 



(57) 



with the counting parameters A = {Xei ^n}, the fluctuating instantaneous currents j(i) = {j£(i), jAr(t)}, and where 
(•)a denotes the statistical average in the stationary state of affinities A. This generating function can be obtained 
as solution of the eigenvalue problem: 

^a* = -Qa(A)* (58) 
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for the modified master-equation operator: 

L x =c- XG Lc XG (59) 

where G = J Q j(r) dr are observables called Helfand moments in transport theory and such that their time derivative 
is equal to the instantaneous current j(t) — dG/dt [35|. We have the following: 



Theorem. The modified operator Eq. H59\) corresponding to the master equation H28\) obeys the symmetry relation: 

I) - 1 eL A (e !) $) = ii_ A $ (60) 

where <I>(N) denotes some arbitrary function, O is the involution performing the time-reversal transformation: 

0$({N a }) = $({7V qT }) so that 9 2 = 1 (61) 
and r)(N) is the multiple Poisson distribution: 

tj(N) =l[ c -( N ^^l— with (N a ) = f R (v a ) Ar 3 Av 3 (62) 
l\ a . 

a 

given in terms of the Maxwell- Boltzmann distribution /r(v) of some reference reservoir, for instance the right-hand 
reservoir at the temperature Tr and density 

B. Proof of the theorem 

1. Open system in contact with two reservoirs 

This theorem is established as follows. The first step is to notice that only the flow part of the master-equation 
operator is modified by Eq. (l59l) : 

L\ E ,\ N = L™, A N + L ( 63 ) 

Indeed, the operator L c due to the binary collisions acts locally exchanging particles among all the cells corresponding 
to the same position. Since the observables (|55[) and (|56[) are non zero for particle exchanges at the boundary with 
the left-hand reservoir, there is no modification of the binary-collision operator. 
Now, the binary-collision operator always obeys the symmetry relation: 

^'ei (e?)$) = L ct $ (64) 

as shown in Appendix [S] by using the property that the lowering and rising operators (|25l) act as follows on the 
product of the multiple Poisson distribution (|62[) with an arbitrary function $: 

E- 1 r 1 ® = 1 ^-r ] E- 1 $ (65) 

^^^^Ti 11 ^ 1 * (66) 

On the other hand, the flow operator (l4"2"j) is modified by the transformation ([55)1 because it contributes to the 
exchange of particles with the left-hand reservoir where the currents (|55|) and (|5"6"|) are defined. For particles at the 
energy e = mv 2 /2, the flow operator (1421) is modified into 



Ar 



i-i 



(e" A - l){N+} h $ + Y,(E+lE- + \, + - l)JV i)+ # + (E+X - 1)^,+$ 



i-i 



(c+ x E+X - l)Ni,-$ + Y,^tXi,-K- - + i E 7X - 1 )( iV ->R $ 



(67) 
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with the counting parameter A = eA^ + Ajv- As aforementioned, the multiple Poisson distribution (|62[) is symmetric 
under time reversal, Or/ — 77, where the time-reversal transformation acts as QNi ± = exchanging the velocities 

±v at every position i = 1,2, /. Moreover, we have the properties (|27p . (|65l> . and (|66|) for the lowering and rising 
operators. Therefore, the transformation in the left-hand side of Eq. (|oT)]) has the effect of attributing the terms 
of positive velocity to those of negative velocity and vice versa. Furthermore, the boundary term with e~ x (N + )l is 
transformed into e~ A (iV+)L/ (Ni t -)-R and the boundary term with e +A into e +A (A^i,+)R,. Now, the left-hand reservoir is 
at the temperature Xl and density riL while the right-hand reservoir is used as reference so that (Ni t ±)n = f^Ar 3 Av 3 
for all i = 1,2,...,/ according to Eq. (|62p. whereupon we have the equalities: 

c- A <g±>± = c^ A (68) 

c +X (N h+ ) R = e-^- A > (N + ) L (69) 

with the affinity A defined by Eq. (|43|) . Consequently, the two boundary terms with the counting parameter A may 
exchange their role if we carry out the transformation A — >• A — A. We thus obtain the symmetry relation of the flow 
operator: 

Combining the results and (|701) . the symmetry relation ([SO]) is finally proved. 

2. Open system also in contact with a thermahzing surface 

Here, we consider the case where the system is open and includes a thermahzing surface described by the operator 
(1451) . Since energy exchange may occur at the surface, we have to introduce an extra counting parameter X% cor- 
responding to an instantaneous current similar to Eq. (|55[) but defined at the surface. Accordingly, the operator is 
modified to 



£ A|* = E 5 -.o E E Av " ^ P^W) [e A|(£ '- £) 4 + v' E- 1 l] N rv $ (71) 

where e = mv 2 /2 and e' — mv' 2 /2. The identity 

V-'OLls (9 7?$) = L% s $ (72) 
is obtained provided that the gas-surface interaction kernel satisfies: 



\v z \p(-v'\ - v) 



e lAS-PRjl«-« J (73-) 



This is always guaranteed thanks to the property of reciprocity (1201) if the affinity associated with the thermahzing 
surface is defined by 

A% = (3 K - p s (74) 

where (3 S — (k-sTs)^ 1 characterizes the Maxwell-Boltzmann distribution (fT2|) at the surface temperature T s . 

For an open system in contact with two reservoirs and a thermahzing surface, we have thus proved the symmetry 
relation ([50]) with the counting parameters A = {A#, Aat, A|;} and the affinities A = {Ae 1 An, A^} given by Eqs. (|13p . 
P]>, and (j?3]>. 

Such symmetry relations extend from the coarse-grained master equation to the fluctuating Boltzmann equation as 
long as this latter is the diffusive limit of the former together with their respective modified evolution operators. 

C. The current fluctuation theorem 

The consequence of the theorem is that the real eigenvalues of Eq. (|58[) and, in particular, the leading eigenvalue 
giving the generating function (1571) also have the symmetry of the modified operator so that we get the time-reversal 
symmetry relation for the cumulant generating function: 

Qa(A) = Qa(A — A) (75) 
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If we denote by Pa(J,£) the time-dependent probability that the instantaneous currents averaged over a time 
interval t take the values 



1 « 



J = ~ / j(r)dT (76) 



the fluctuation theorem for the currents is obtained: 



P A (-J,i) 



exp (A • J t) (77) 



in the long-time limit for the fluctuating Boltzmann equation. The equivalence between Eqs. (|75l) and (|77l) is demon- 
strated by noting that the cumulant generating function can be expressed in terms of the probability of the fluctuating 
currents as 

Q A (A) = lim --\n [ P A (J,t) exp(-A ■ Jt) dJ (78) 

t— >oo t J 

so that Eq. (|77|> yields Eq. (|75jl and vice versa. 

As a consequence of the current fluctuation theorem, the generalizations of Green-Kubo formulae and Onsager 
reciprocity relations obtained in Ref. [19] for the nonlinear response coefficients extend to the nonequilibrium flows of 
dilute and rarefied gases ruled by the fluctuating Boltzmann equation. 



VI. CONCLUSIONS 



In this paper, a time-reversal symmetry relation has been established for all the random currents of particle and 
energy carried by the particles of dilute or rarefied gases flowing in open systems ruled by the fluctuating Boltzmann 
equation. For this purpose, the fluctuating Boltzmann equation is obtained as the Langevin equation corresponding 
to a coarse-grained master equation for the random numbers of particles in cells partitioning the one-particle phase 
space, in the limit of arbitrarily small cells. The coarse-grained master equation has been obtained not only in the 
bulk of the fluid but also at its boundaries where the nonequilibrium constraints are imposed. These constraints 
are characterized in terms of the thermodynamic forces or affinities defined from the differences of temperatures or 
particle densities between gas reservoirs or, possibly, thermalizing surfaces. 

The time-evolution operator of the coarse-grained master equation is composed of several terms corresponding to the 
different processes taking place in the bulk or at the boundaries. This linear operator is modified by introducing the 
observables counting the transfers of particles and energy at the boundaries and the associated counting parameters. 
The cumulant generating function of the fluctuating currents is given as the leading eigenvalue of this modified 
operator. The central result is obtained by showing that the modified operator obeys the symmetry relation (|6TJ)) 
involving the counting parameters and the affinities. As a consequence, the cumulant generating function has the 
symmetry (1751) and the fluctuation theorem ([77)1 is thus established for all the currents. These relations find their 
origin in the time-reversal symmetry of the underlying microscopic dynamics. 

The fluctuation relation (|77l) implies the non-negativity of the entropy production defined as the sum of the affinities 
multiplied by the average values of the currents flowing across the open system. Furthermore, the fluctuation relation 
also implies the Green-Kubo formulae and Onsager reciprocity relations as well as their generalizations beyond linear 
response theory [l9j . which thus hold for dilute or rarefied gases ruled by the fluctuating Boltzmann equation, as the 
consequence of the symmetry relations here established. 

These results concern in particular the numerical simulation of dilute or rarefied gases by methods such as the 
Direct Simulation Monte Carlo method if the enunciated assumptions and, especially, the properties of reciprocity 
(fTUj) and (H2) are satisfied 

At equilibrium where the affinities vanish, the fluctuation relation (I77|) reduces to the expression of the principle of 
detailed balancing, according to which the probabilities of opposite fluctuations are equal. Out of equilibrium, a bias 
appears between these probabilities. Since the ratio of the probabilities of opposite fluctuations goes exponentially 
with time and the affinites, one of both probabilities may soon become so tiny that it would appear negligible, as 
nonequilibrium constraints increase. This is the case in regimes very far from equilibrium where the fluxes tend to 
become unidirectional. In extreme regimes for instance with shock waves, the fluctuation relation thus expresses the 
quasi full irreversibility of the flow due to the extreme rarity of fluctuations opposite to the average fluxes. It is only 
relatively close to equilibrium that the opposite fluctuations are not negligible. 

In the present paper, the results have been obtained for gases composed of a single atomic species. However, the 
results extend to polyatomic molecular species, gas mixtures, as well as reacting dilute or rarefied gases, as long as 
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the property of microreversibility may be assumed [36|, |37[ . The scope opened by these results is very broad since 
we may also envisage their extension to gases ruled by the relativistic Boltzmann equation or to gases of bosons or 
fermions ruled by the quantum Boltzmann equations [H, H§| . 

Acknowledgments 

The main result of this paper has been presented at the conference "Boltzmann equation: mathematics, modeling 
and simulations" in memory of Carlo Cercignani at the Henri Poincare Institute, Paris, February 9-11, 2011. This 
research has been financially supported by the Belgian Federal Government (IAP project "NOSY"). 

Appendix A: Time-reserval symmetry of the binary-collision operator 

In this appendix, the time-reversal symmetry relation (|64[) is proved for the binary-collision operator 



L c 4> = E W \»p° (E^ 1 E- 1 E+ 1 E+ 1 - l) N p N a $ (Al) 

Using Eq. (f27|) . its adjoint is given by 

L°t$ = J2 W^ pa N p N a ( k E; 1 E- 1 E+ 1 E+ 1 - l) $ (A2) 

Xppcr 

Using the time-reversal operator (|rJTj) and Eqs. (|65p - (jf)6]) , the expression in the left-hand side of Eq. is transformed 
as follows: 

r^G L c (6 r) *) = r 1 E W ^p° E$E-iE$E£ A> Y ffT r, $ - E W? ppa A> NaT <j> 

Xfipcr Xfipcr 
Xppa " Xfipcr 

Now, summing over all the cells a gives the same result as summing over the time- reversed cells a /T , so that we get 
rj-'S L c (6 ry $) - E Kw** W^^S Np N ° V^^A ^ $ ~ E ^pW^t N x N p $ (A4) 



By using the property of reciprocity (|41j) with respect to the equilibrium state of the right-hand reservoir, the first 
term becomes 

^9^(8^)=^ W^ pa N p N a E- 1 E- 1 E+ 1 E+ 1 * - £ w£ ffTAT/lT N x N p $ (A5) 

Xfipcr Xp,pcr 

By using Eqs. (|38]) - (l39|) to transform the second term, we finally obtain 

rj-'QL (9^)=E - l) $ (A6) 

A/^pcr 

which is identical to the adjoint operator (|A2[) . Q. E. D. 
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